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ABSTRACT 

We present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z > 1 3CR sources, from the Herschel 
GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we 
perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the 
star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution 
from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG- 
strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black 
hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation- 
based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of 
equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation 
(negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether 
their supermassive black holes are accreting or not. 
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1. Introduction 


The understanding that most (if not all) galaxies in the Uni¬ 
verse host a supermassive black hole (SMBH) is among the 
most important findings of modern astronomy. The growth of 
a SMBH through mass accretion generates large amounts of en¬ 
ergy during a phase in the evolution of the galaxy known as an 
active galactic nucleus (AGN) phase. Although there is a dif¬ 
ference of a factor of ~ 10 9 in their physical size scales, the 
SMB Hs and their host galaxies exhibit strong scaling relations 

fe.g.lMagorrian et alJll998LlTremaine et alJl2002l : lGultekin et all 

20091) . suggesting a link between the growth of the SMBHs and 
that of their host galaxies. Moreove r, both these processes are 
thought to peak at redshifts z ~ 2 (e.g. lHopkins & Beacoml2006l : 

I Alexander et al.l [20081) . The symbiosis of black hole and global 
galaxy growth is intriguing because of the possible feedback ef¬ 
fects: positive (AGN inducing star formation) and/or negative 
(AGN quenching of star formation). These feedback processes 
are of paramo unt importance for our understanding of galaxy 
formation (e.g. ICroton et al .1120061 iHonkins et akll2008l) . How¬ 
ever, neither the feedback mechanisms nor the overall impact of 
feedback on the host galaxies is known. Other big unknowns are 


* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


the duration and frequency of AGN accretion and host galaxy 
star formation phases. 

High redshift radio-loud AGN (P I 4 ghz > 10 25 W Hz ' and 
Z > 1) provide a unique opportunity to probe the interplay 
between the growth of the black hole and the hosting stel¬ 
lar bulge. They are invariably associated with massive galax¬ 
ies having M sfP u ar > 10 11 M n dBest et al.l 1 19981: ISevmour et al.l 
120071:iDe Breuck et al.12010), and have edge-brighte ned, double- 
lobed, FRII morphologies ( Fanaroff & Rilevl 119741 ) that permit 
estimates of the duration of the episode of strong AGN activity. 
In addition to being used in studies of massive galaxy evolu¬ 
tion, radio-loud AGN are being used extensively in unification 
studies, where, distant radio-loud galaxies and quasars are be- 
lieved to make up one and the same popu lation ( Ba rthellll989b 
lAntonuccil 119931: lUrrv & Padovanil 1 19951) . hence have equally 
massive hosts. Ultraviolet or visible opaque circumnuclear dust 
is an essential element of this scenario; distant 3CR quasars and 
radio galaxies are indeed luminous mid-infrared (MIR) emit- 
ters ( Siebenmorgen et akl 120041: lHaas et akl 120081: iLeinski et al.1 
120101 ) . 


It has long been suspected that hosts of powerful high- 
redshift radio-loud AGN un dergo episodes of vig orous (dust ob¬ 
scured) star formation (e.g. [Archibald et al.llToOli) . Huge reser¬ 
voirs of molecular gas, have been deduced in several objects 
from submillimetre (submm) studies (e.g. [Reuland et al. 200 4). 
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Such studies were mainly based on one submm flux measure¬ 
ment and were limited to the highest redshift objects for which 
the obscured newborn star radiation re-emitted by the ubiquitous 
cold (30-50 K) dust is redshifted to submm wavelengths. How¬ 
ever, quantification of the cold dust emission (e.g. constraining 
the cold dust temperature) requires sampling the full rest-frame 
infrared-submm spectral energy distribution (SED) of the stud¬ 
ied objects. Earlier far-infrared (FIR) studies failed to provide 
strong constraints on the FIR properties for relatively large sam¬ 
ples of powerful radio-loud AGN because of their small detec- 


age ( 

Heckman et al. 1992^ 

Hes et al. 1995; Meisenheimer et al. 

2001 

; Siebenmorgen et al. 

2004; Haas et al. 2004; Clearv et al. 

2007 

). 


The Herschel Space Observatory dPilbratt et al.ll2010h . with 
its unprecedented FIR sensitivity and wavelength coverage, ex¬ 
plored terra incognita (caelum incognitum...) allowing studies 
which have revolutionized the understanding of the connec¬ 
tion between AGN and star formation activity. Several stud¬ 
ies utilizing deep X-ray and Herschel data revealed that the 
hosts of moderately luminous radio-quiet AGN out to z ~ 3 
form star s at rates comparable to the general non-AGN pop¬ 
ulatio n (IShao et alJ 1201 Ot iMullanev et al.1 1201 2l iRosario et al.1 
12012). For high (radio- quiet) AGN luminosities (L x > 10 44 erg 


s 1 ), Page et af] d2012l) reported suppression of star formation, 
consistent with the expectations from theoretical models, while 
lHarrison et al.l (i2012l) found no clear evidence of suppression of 
star formation by extending the analyses to samples larger by an 
order of magnitude. Moreover, at the highest AGN luminosities 
(in excess of 10 46 erg s” 1 ), recent studies based on decompo¬ 
sition of the IR emission to AGN and star formation contribu¬ 
tions, have shown star formation rates (SFRs) of the order of 
several hundred solar masses per year in the hosts of some of the 
most powerful radio galax ies (Barthel et al. 20121 - B12 here¬ 
after; ISevmour et a 1 .1 201 2t Drouart et al .1120141) and (radio-quiet) 
quasars (Feinsk iet al.ll2013T 20141) . 

In order to quantify the energetics of AGN at the peak of 
their activity as well as their star formation characteristics, we 
obtained five-band Herschel photometry of the 3CR sample us¬ 
ing the Photodetector Array Camera (PACS) at 70 and 160 pm 
and the Spectral and Photometric Imaging Receiver (SPIRE) at 
250, 350, and 500 pm on-board the Herschel Space Observatory. 
The first results, dealing with 3 archetypal objects of that sample 
were presented in B12. Here we analyse the FIR properties of the 
complete (flux-limited) sample of objects spanning the redshift 
range 1 < z < 2.5. This paper is organized as follows. Section[2] 
describes the sample selection, the data obtained, and the steps 
used for measuring flux densities in the five Herschel bands. Sec¬ 
tion [3] addresses the procedure for fitting the observed IR SEDs 
of the objects. Results are then presented and discussed in Sect. [4] 
and Sect. 0 respectively, and the paper is briefly summarized in 
Sect. [6] Throughout this paper we use a flat cosmology with H<> 
= 70 km s" 1 Mpc -1 and fl,\ = 0.7, an d we follow the conversion 
in lKennicutd (fl998) (which assumes a Salpeter 1955 initial mass 
function) when deriving SFRs. 


2. Data 

2.1. Sample selection 

With this study, we target the well known, complete flux-limited 
sample of the brightest (F 1 7 g MHz >10 Jy), high-redshift (z > 1) 
radio-loud AGN sample in the northern hemisphere: the Revised 
Third Cambridge Catalogue of radio sources (hereafter 3CR; 



Fig. 1 . Observed radio (178 MHz) luminosity as a function of redshift 
for the 3CR sample considered in this work. Circles (red) indicate radio 
galaxies, and squares (blue) indicate quasars. The plus symbols mark 
the objects detected in at least three Herschel bands (typically the two 
PACS and the SPIRE 250 pm bands. 


ISninrad et al.lll985l) . The extremely high luminosities (Fig. [TJ of 
these double-lobed radio galaxies (RGs) and quasars (QSRs) are 
produced by some of the most powerful accreting SMBHs. The 
low-frequency (178 MHz) radio selection ensures no bias with 
respect to orientation: the steep-spectrum lobes of radio-loud 
AGN emit optically thin and isotropic synchrotron radiation, 
making the 3CR sample ideal for testing the orientation-based 
unification scenario of radio-loud AGN. As shown in Fig.Q] both 
the RGs and QSRs are homogeneously distributed in redshift. 
The z > 1 3CR objects show mostly FRII morphologies; this is 
well-established from, for instance, high-resolution VFA maps. 
Compact, presumably young, morphologies dO’Dea|[l998l) are 
also found within the sample. 

The z > 1 3CR sample is (spectroscopically) completely 
identified using 3 to 5 m-class telescopes in the 1960s-1980s 
(iSpinrad et al.lll985l) . The objects in the sample almost univer- 
sally accrete at high E ddington rate, i.e. in quasar-mode (e.g. 
iBest & Heckmanl2012l) . The total number of objects in the z > 1 
3CR sample is 6dfl The highest redshift 3CR source is 3C 257 
(z. = 2.47). Two z > 1 3CR sources, 3C 287 and 3C 300.1, 
have been observed in other Herschel observing modes, and thus 
have been dropped from this work. The remaining 62 sources, 37 
RGs and 25 QSRs, comprise the Herschel sample studied in this 
work. An overview of selected properties is provided in Tahiti] 

The high-z 3CR sample has been observed with many spac e 
telescopes (including Hubble , Spitzer, and Chandra ; IBest et al.l 
[l998t lHaas et al. 1120081 [Ldnski et al.ll20l(it IWilkes et alJl2013l) : 

objects from this sample represent landmarks in our study of ac¬ 
tive galaxies through cosmic time. 

2.2. Herschel photometry 

The data for this work were obtained as part of our Herschel 
Guaranteed Time project The Herschel Legacy of distant radio- 
loud AGN (PI: Barthel, 38 hours of observations). Five objects 
(see Table Q} were observed as part of another Herschel pro¬ 
gramme (PI: Seymour). The raw data for these objects were re- 

1 The sample includes two 4C objects, 4C 13.66 and 4C 16.49, which 
formally match the selection criteria of, and are included in, the 3CR 
sample. 
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trieved from the Herschel Science Archive (HSA), and the data 
reduction was performed as detailed below. 


2.2.1. PACS 

Photometric observations were carried out with PACS 
JPoglitsch et al.l feoiot) in the scan-map observational mode, 
both in the blue (70 pm, 5" angular resolution) and in the red 
(160 pm, 11" angular resolution) bands. A concatenated pair 
of coextensive scan maps at two different orientations was 
obtained for each source. Data reduction was performed within 
the Herschel Interactive Processing Environment (HIPE, lOtt) 
12010L version 11.0.0), following the standard procedures for 
deep field observations. Maps were created by employing the 
high-pass filtering method, using an appropriate source masking 
step to avoid significant flux losses due to the high-pass filter. A 
first data reduction resulted in a preliminary map, created after 
combining the individually (for each scan orientation) processed 
scan maps. Source masking was performed by hand, using the 
preliminary created map as an input. This method in particular 
allowed us to minimize the flux losses of the observed sources 
dPonesso et al.il20f2h . The final data reduction and mosaicking 
were then performed using the mask generated in the previous 
step. 

Photometry (using appropriate aperture corrections) was per¬ 
formed within HIPE, using the annularSkyAperturePhotometry 
task. Apertures of 6" and 10" radius for PACS blue and PACS 
red, respectively, were centred on the known radio core posi¬ 
tion of the object in the map. PACS maps suffer from correlated 
noise, thus pixel-to-pixel variations cannot yield robust photo¬ 
metric uncertainties. Instead, we opted for the well-established 
procedure of placing apertures at ran dom positio ns on the sky 
( Lutz et af] 1201 it [Popesso et al.ll2012l) . Following iLeinski et al.l 
( 2013 ). we placed 500 apertures of 6" (for blue) and 10" (for 
red) radii at locations avoiding the noisy edges of the map, re¬ 
quiring that the central pixel of the random aperture has at least 
75% of the integration time of that of the source of interest. The 
resulting distribution of the flux densities measured in these 500 
apertures was then fitted with a Gaussian, and the sigma value 
of the Gaussian was taken to be the lcr photometric uncertainty 
of the map. Measured flux densities and associated lcr uncertain¬ 
ties, together with 3cr upper limits for the non-detections are pro¬ 
vided in Table [2] The PACS photometric uncertainties provided 
in Table[2]do not include the 5% uncertainty on the absolute flux 
calibration dBalog et al.ll2013 j). Postage stamps of the resulting 
PACS maps, centred on the radio position of the objects, are in¬ 
cluded in Appendix[E] 


2.2.2. SPIRE 

SPIRE (iGriflfin et al.l l2010l ) photometric observations were car¬ 
ried out in small scan-map observational mode, at 250 (18.2" 
angular resolution), 350 (24.9" angular resolution) and 500 pm 
(36.3" angular resolution). Data reduction was performed in 
HIPE following standard procedures for SPIRE data. Source 
extraction on the fully reduced m ap was performed usin g the 
sourceExtractorSussextractor task rfSavage & Qliveifc007t) . Ex¬ 
tracted sources located within half the Full Width at Half Max¬ 
imum (FWHM) of the given SPIRE array (measured from the 
radio core position of the sources) were selected as tentative de¬ 
tections. 

While SPIRE does not suffer from correlated noise, SPIRE 
observations are dominated by confusion noise, of the order of 


6-7 mJy beam -1 , as e stim ated from deep extragalactic observa¬ 
tions dNguven et al.ll2010l) . Our adopted procedure for the de¬ 
termination of the photometric uncertainties in the SPIRE maps 
is fully described bv ILeinsk i et al.l (1 201 3l). whi ch i n turn follows 
the pro cedures presented bv lElbaz et al.l (1201 ll) and lPascale et al.l 
(1201 ll) . Initially, an artificial source free map was created by re¬ 
moving the extracted sources in a SPIRE map from the map it¬ 
self. Then, the pixel-to-pixel rms in this source free map was 
calculated using a box centred on the nominal position of the 
target object. The size of the box was chosen as a compromise 
between avoiding the noisy edges of the SPIRE map and obtain¬ 
ing proper statistics of the immediate environment of the target 
object. 

As indicated in Table [2] three SPIRE detections are for¬ 
mally below the estimated 3cr values. These particular measure¬ 
ments were included in the subsequent analyses because, upon 
visual inspection, they showed obvious emission at the known 
position of the target object. The availability of ancillary multi¬ 
wavelength data at shorter wavelengths (to check for source con¬ 
fusion), and the understanding of the overall shape of the ob¬ 
ject’s SED, further support the inclusion of these flux densities 
in the subsequent analyses. While the formal signal-to-noise ra¬ 
tio of one of these three detections is very close to three, the 
other two detections do not reach this ratio only because the as¬ 
sociated SPIRE maps are less clean than other maps in the sam¬ 
ple, leading to significantly larger photometric uncertainties. The 
SPIRE 500 pm photometry should be considered tentative be¬ 
cause the beam at this particular wavelength is large, and unde¬ 
tected sources in the region surrounding the AGN may contribute 
to the measured flux density. The SPIRE photometric uncertain¬ 
ties provided in Table [2]do not include the 4% uncertainty on the 
absolute flux calibration (iBendo et al.ll2013l ). Postage stamps of 
the resulting SPIRE maps, centred on the radio position of the 
objects, are included in Appendix[E] 


2.3. Supplementary Data 

The FIR photometry of all 3CR sources in our work was 
supplemented with MIR photo metry obtained with the Spitzer 
Space Telescope (I Werner et al.l 120041) during three Spitzer GT 
observing programmes (PI: G. Fazio) in six bands, using the 
instruments IRAC dFazio et alJ I20041). IRS-16 pe ak-up array 
dHouck et al.1 120041) . and MIPS dRieke et~akl 120041) . Details on 
the Spitzer data reduction and photometry have previously been 
published by lHaas et all d2008l) . Table [2 lists the Spitzer pho¬ 
tometry. When available, additional 850 pm data were collected 
from the literature. The 850 pm emission in quasars can be heav¬ 
ily contaminated by synchrotron contribution (see Sect. 13.31) . The 
quasar 850 pm thermal flux densities utilized in this work were 
taken from lHaas et all (| 2006|). To obtain the 850 pm thermal 
flux densities of quasars, Haas et all J2006) extrapolated the syn¬ 


chrotron contribution at 850 pm using the measured radio core 
flux densities, and subtracted it from the total flux density at 
850 pm. Table[3]lists the radio galaxy and quasar thermal submm 
flux densities used in this work. 


3. Spectral energy distributions 

3.1. Fitting components 

The estimation of physical properties for the active galaxies was 
performed using an SED fitting technique. Our fitting routine is 
based on a combination of several distinct components, responsi¬ 
ble for the emission from active galaxies in different wavelength 
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Fig. 2. IR spectral energy distributions (SEDs, solid black) for two representative objects from this work. Open circles show the photometric data. 
Error bars correspond to 1 <x photometric uncertainties. Arrows indicate 3 tr upper limits. Left panel: 3C 454.1, a radio galaxy at z=1.84. The three 
components used to fit the SEDs of radio galaxies account for emission from host galaxy (old) stars (dash-dotted red), from an AGN-heated torus 
(dotted green), and from dust heated by star formation (dashed yellow). Right panel: 3C 205, a quasar at z=1.53. The multi-component approach 
used to fit quasars accounts for emission from hot (1300 K) dust (long-dashed blue), from the AGN-heated torus (dotted green) and from the star 
formation heated dust (dashed yellow). 


regimes. Below, we describe this multi-component approach of 
fitting the observed IR SEDs. 

The presence of circumnuclear dust surrounding the broad 
line regions in AGN and blocking their UV/visible emission 
is central to orientation-bas ed unification of powerful FRI I 
radio galaxies and quasars (Barthel 1989; Antonu ccil 1 19931) . 
Given its proximity to the AGN, the emission from this AGN- 
illuminated war m dust peaks at rest-fram e MIR (10-40 pm) 
wavelengths (e.g. iRowan- Robinsonl[l995 ). Spitzer photometric 
and spectroscopic data have shown that the majority of high- 
z 3CR objec ts are luminous MIR emitters (lHaas et al.l 2008; 
Leins ki et al.1 [2010), with observed MIR luminosities L 1 5;jrn 
much higher than 8 x 10 43 ergs s 1 , the value separating hid¬ 
den quasars from mid-IR weak radio galaxies at intermediate 
redshifts dQgle et al.ll2006l) . There exists broad agreement that 
the AGN-heated nuclear dust is mainly located in clumps which 


AGN torus (e.g. Nenkova et alJ 200^ 

Kuraszkiewicz et al. 20031 

Honig et al. 2006; Schartmann et alJ 

2008). To account for the 


emission from the AG N heated dust, we cho se the library of 
torus models from lHonig & Kishimotol d2010t) . The parameters 
considered when creating this library include the (1) radial dust 
distribution of dust clumps; (2) geometric thickness of the torus; 
(3) number of clumps along an equatorial line of sight; (4) op¬ 
tical depth of the individual clumps; and (5) outer radius of the 
torus. There are 240 sets of parameters in the library, each com¬ 
puted for seven viewing-angles ranging from face-on (i = 0 ) to 
edge-on (i = 90°) in steps of 15 degrees, resulting in a total of 
1680 torus models. Detailed information on the model parame¬ 
ters and t he adopted strategy in gen erating the tori SEDs is pro¬ 
vided bv lHonig & Kishimotol ( 2010 ). In addition to the parame¬ 
ters listed above, the overall flux normalization of the model is 
another free parameter throughout the fitting procedure (outlined 
below). 

The rest-frame FIR emission (40-500 pm) is largely gener¬ 
ated by cold du st, heated by star formation on kpc-scales in the 
AGN host (e.g. Rowan-Robinsonlf 1 9951; ISchweitzer et alJl2006t 


iNetzer et al.l 12007 ). Following these authors, we interpret any 


FIR emission (in excess of the emission of the AGN-heated dust) 
as being powered by star formation, and we represent it with an 
optically thin modified blackbody component, i.e., a blackbody 
modified by frequency-dependent emissivity, given by 

S v ocB v (T)\P. (1) 

We reduce the number of free parameters in this component by 
fixing the dust emissivity index f5 to a typical value of 1.6 as 
found in studies of high-;; AGN (e.g. lBeelen et al,f2006l) . The re¬ 
maining two free parameters here are the cold dust temperature 
and the flux normalization of the modified blackbody compo¬ 
nent. The use of a m odified blackbody, as opposed to starburst 
templates (e.g. lDrouart et al.l2014l) . might slightly underestimate 
the star formation luminosities because one misses the starburst 
MIR emission, but offers the unique possibility of constraining 
the cold dust temperatures. 

The two components described above feature in the SED fit¬ 
ting of both radio galaxies and quasars. We include additional 
SED components to the fitting depending on the type of the stud¬ 
ied object. For radio galaxies, we added a blackbody component 
peaking in the near-IR to account for the emission from the old 
stellar population in the AGN host (e.g. ISevmour et al.ll2007h . 
The temperature of the blackbody and its flux normalization are 
the two free parameters for this SED component. For quasars, we 
added a blackbody component to account for the hot (graphite) 
dust close to the sublimation temperature. This component is 
often empirically required to fit the observed SEDs of quasars 
(e.g. iMor & Netzerll2012l; ILeipski et~aDl2013l) . Following these 
authors, we fixed the temperature of the blackbody to 1300 K, 
leaving its flux normalization as the only free parameter during 
the fitting. Such a component is also needed in the fitting of the 
SEDs of a few radio galaxies (see also B12) whose observed 
photometry in the NIR/MIR could not be well represented with 
the components described above. These radio galaxies, indicated 
in Table 01 might be viewed along lines of sight at which the nu¬ 
clear region is only partly obscured, thus resulting in somewhat 
elevated MIR luminosities. The inclusion of the hot dust com¬ 
ponent to the SEDs of some of the radio galaxies might lower 
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the estimates of the mass of the evolved stellar populations, but 
this is outside the scope of this work. For the quasars we also 
considered an additional power-law component representing the 
emission from an accretion disk in the UV/visible. However, as 
demonstrated in Appendix [A] the inclusion of this power-law 
component had little influence on the results obtained from the 
FIR part of the SED, therefore it was excluded from the fitting 
procedure. 


3.2. Fitting procedure 


While our physically motivated fitting approach results in a close 
approximation of the observed SEDs of the sample objects, it is 
not primarily designed to yield precise models of their SEDs. 
In particular, we are not interested in constraining the proper¬ 
ties of the dusty torus with our multi-wavelength broad-band 
photometry. This kin d of analysis remains challenging even at 
lower redshifts (e.g. iRamos Almeida et al. 2009 ). Torus mod¬ 
els are used to separate AGN-heated dust emission (peaking 
in MIR) from star-formation-heated dust emission (peaking in 
FIR), and to determine for the first time the star-formation- 
dominated FIR energetics of the high-’ 3CR sources. When fit¬ 
ting the observed SEDs of our objects, we used a chi-square 
minimization technique based on the MPFIT routine (Markwardtl 
2009) . In practice, we started with a torus model from the library 
of Honig & Kishimoto (2010) and added a linear combination of 
the remaining SED components (depending on the object type) 
to minimize the overall chi-square. We repeated the procedure 
for each torus model in the library. 


Example best-fit SEDs, along with their individual SED 
components, are shown in Fig. [2] At the redshifts of our sam¬ 
ple, the PACS 70 pm band is crucial in our adopted fitting ap¬ 
proach, as it strongly constrains the longer wavelengths of the 
torus emission. On the other hand, the SPIRE 250 pm band is 
the most important measurement for constraining the compo¬ 
nent representing the cold dust emission. Therefore, our fitting 
approach was applied to all objects that are detected in at least 
three Herschel photometric bands (typically the three shortest 
Herschel bands). These objects are homogeneously distributed 
in the redshift range (1 < z < 2.5) studied in this work (Fig. [TJ- 
Best-fit SEDs for these objects, together with images centred on 
the radio positions of the AGN, are presented in Appendix [C] 
and in Appendix[E] respectively. Occasional SED mismatches at 
observed-frame 16 pm and/or 24 pm are most probably due to 
luminous polycyclic aro matic hydrocarbon emission and/or th e 
10 pm silicate absorption (lHaas et al.ll2008ULeipski et al.l l2010). 
At FIR wavelengths, the fixed f3 approach may be the main rea¬ 
son behind the failure of the fitting routine to exactly reproduce 
some SPIRE data points. 


The physical parameters constrained by our fitting method 
include IR star formation and IR AGN luminosity, and the tem¬ 
perature and mass of the cold dust component (Fig. [3] Table |4j. 
Uncertainties in the derived parameters were calculated by re¬ 
sampling the observed SEDs, allowing the individual photo¬ 
metric measurements to vary within their lcr ranges of uncer¬ 
tainty. More precisely, we generated 100 mock observed SEDs 
and studied the distributions of the parameters derived from the 
corresponding best fits. We inspected the best fits to the mock 
SEDs to confirm their overall quality in the different wavelength 
regimes. From the distributions, we retained the median values 
as the best estimates of the parameters, and the 16th-84th per¬ 
centile ranges as their associated uncertainties (which in case of 
a Gaussian distribution would correspond to ±lcr values). 


For objects with fewer than three Herschel detections, we 
estimated upper limits for the IR star formation and IR AGN lu¬ 
minosities using two different approaches. In the first approach, 
we fitted the (longer wavelength) Herschel upper limits using 
a modified blackbody with fixed f3 - 1.6 value and fixed cold 
dust temperature (7jj ust = 37 K), typically found from the fits of 
the objects detected in at least three Herschel bands (see below). 
We then integrated under the blackbody component to estimate 
upper limits for the star formation luminosities. In the second 
approach, we took the 70/160/250 pm upper limits as tentative 
detections, and estimated upper limits for the star formation lu¬ 
minosities using the procedure adopted for the objects detected 
in at least three Herschel bands. Both approaches yielded simi¬ 
lar results (within 10%) for the star formation luminosities, but 
we retained the second approach because it allowed us to es¬ 
timate the IR AGN luminosities for objects with only PACS 
70 pm detections. Best-fit SEDs for objects with fewer than 
three Herschel detections, together with images centred on the 
radio positions of the AGN, are presented in Appendix[D]and in 
Appendix[E] respectively. Systematically demanding three Her¬ 
schel detections when fitting the observed photometry regardless 
of redshift means that the FIR results for the objects detected in 
only the two PACS bands are treated as upper limits. However, 
depending on the object’s redshift, the PACS 160 pm band alone 
may probe the peak of the cold dust emission, allowing robust 
constraints for the physical parameters estimated in this work. 


3.3. Synchrotron contribution 


Earlier submillimetre and millimetre studies of high-redshift 
3CR sources presented clear evidence for a synchrotron contri- 


bution to the observed flux densities (van Bemmel et al 


lArchibald et alJ 1200 it IWillott et al.1 120021: Haas et al.1 


1998; 

20061 ) . 


While negligible for radio galaxies, extrapolation from core ra¬ 
dio data shows that synchrotron emission can account for up 
to 8 0% o f the observed submm flux densities from quasars 
(lHaas et al.ll2006f) . However, the longest rest-frame wavelengths 
probed by the SPIRE 500 pm band in our study are ~ 200 pm: 
emission at these wavelengths is completely dominated by dust 
and therefore free from any synchrotron contribution. The only 
exception is 3C 418. This source appears to be flat-spectrum- 
core-dominated (in an otherwise steep-spectrum selected sam¬ 
ple), and as such its FIR SED is clearly dominated by non- 
thermal (synchrotron) radiation from its core. The power-law 
like IR SED of this source is shown in Appendix IQ 3C 418 was 
removed from the subsequent analyses. 


4. Results 

4.1. Detection statistics 

The Herschel detection rate throughout our sample ranges from 
67% in the PACS 70 pm band to 13% in the SPIRE 500/mi band. 
In particular, 7 objects have robust detections in all five Herschel 
bands. Furthermore, a total of 24 objects are detected in at least 
three Herschel bands, most importantly in the SPIRE 250 pm 
band, which for the highest redshift of our sample (3C 257: z = 
2.47) corresponds to ~ 70 pm rest-frame emissiorQ. Excluding 
3C 418 from these 24 objects, results in 13 radio galaxies and 10 


2 For the median redshift of our sample, (z me d- 1.38), the SPIRE 
250 pm band samples the peak (~ 100 pm) of the typical cold dust SED. 
allowing strong constraints on the modified blackbody component used 
in the fitting approach. 
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Table 3. Objects with significant thermal submillimetre flux densities. 


Object 

Thermal F 85(Vm 
(mJy) 

Reference 

3C 191 

2.95 

1,2 

3C 257 

5.40 

3 

3C 280.1 

2.48 

1,2 

3C 298 

7.25 

1,2 

3C 368 

4.08 

3 

3C 432 

6.33 

1,2 

3C 470 

5.64 

3 

4C 13.66 

3.53 

3 


Re ferences. (DlWillott et alJ(l2002h : (2) lHaas et al](l200fih : 
(3t I Archibald et all <200 111 


quasars. The Spitzer detection rate throughout our sample ranges 
from 94% in the IRS 16 pm and MIPS 24 pm bands to 100% in 
the IRAC 5.8 pm band. Comments on selected individual objects 
are included in Appendix iBl 

4.2. Physical properties obtained from the SED fitting 

Studies investigating the cold dust temperatures in high-! ob¬ 
jects prior to Herschel were often uncertain because they relied 
heavily on observations in a single pho tometric broadband (e.g. 
Benf ord et al.l 11999b iBeelen et al.l 120061) . In our current study, 
cold dust temperatures estimated for the FIR-detected objects 
range from ~ 25 to ~ 45 K (Fig. [3] Table [4}. Radio galaxies 
and quasars span the same range in cold dust temperatures. This 
range is similar to that obtained for z > 5 quasars dLeipski et al.l 
20141) and to that e stimated for distant submm galax ies (e.g. 
Mag nelli et al.l 120121) . Similarly to iLeipski et al.l (1 201 3l) . the in¬ 
clusion of the 1300 K hot dust component in the SED fitting 
lowers the estimates of the cold dust temperatures by -5 K. By 
including the hot dust component, we preferentially select torus 
models that emit more of their energy at longer wavelengths. 
As a consequence, the cold dust components are also shifted to 
colder temperatures. 

Central to the subsequent discussion are the star formation 
luminosities, Lsf, which we computed (or estimated upper lim¬ 
its) by integrating the best-fit modified blackbody components 
from 8 pm through 1000 pm. Figure [3] shows the distribution of 
Lsf for the radio galaxies and quasars detected in at least three 
Herschel bands. Both types of objects show similarly broad dis¬ 
tributions, with many objects having Lsf > 10 12 L e , characteriz¬ 
ing them as ultra-luminous infrared galaxies (ULIRGs) (see also 
Table [4}. The median star formation luminosities of the RGs and 
QSRs plotted in Fig.[3]are 2.0 x 10 12 L 0 and 2.6 x 10 12 L 0 , re¬ 
spectively. Converting the star formation luminosities into SFRs 
using the calibration derived by iKennicutti (1 1998b gives 100 M 0 
yr 1 < SFR < 1000 M 0 yr consistent with SFRs obtained for 
typical sub mm galaxies (SMGs) at comparable redshifts (e.g. 
Ma gnelli et al.ll20l2h . 

To compute the IR luminosities of the components powered 
by the AGN, Lagn, we integrated the best-fit torus component 
for the RGs and the sum of the best-fit hot dust and torus com¬ 
ponent for the QSRf0 between 1 pm and 1000 pm. Figure [3] 

3 A few RGs also require an additional hot dust component to better fit 
their observed photometry. The computation of L AG n for these objects 
takes into account this component as well. 


shows that the RGs and QSRs occupy different ranges in the 
distributions of Lagn, with the distribution for QSRs shifted to 
higher values compared to that for RGs. lHaas et all d2008l) al¬ 
ready established this result by investigating the (energetically 
important) rest 1.6-10 pm wavelength range for the high-! 3CR 
objects, finding the QSRs to be, on average, 3-10 times more 
luminous than RGs. We confirm their finding by including rest- 
frame wavelengths longer than 10 pm. The median IR AGN lu¬ 
minosities of the RGs and QSRs plotted in Fig.|3|are 3.7 x 10 12 
L e and 1.1 x 10 13 L 0 , respectively. 

As outlined in previous sections, we attribute the FIR emis¬ 
sion in excess of the AGN-powered dust emission to emission 
from star-formation-heated dust. Figure |4] shows Lsf as a func¬ 
tion of Lagn for all objects in our sample. The presence/absence 
of correlation between these two parameters depends on both 
redshift and AGN luminosity, and is still debated in the litera¬ 
ture (iLutzl2014l) . Given the data, and taking into account only the 
FIR-detected objects, the two plotted parameters show at most a 
weak correlation, in part introduced by the dependence of both 
Lsf and Lagn on redshift. Moreover, the numerous upper limits, 
together with the fact that both parameters span only a limited 
range (~ 2 orders of magnitude), make it difficult to establish 
any such correlation (or lack of) in our sample. Nevertheless, we 
observe a range of Lsf from weak (if not absent) to very strong, 
coeval with the growth of the black hole. Figure[4]also shows that 
the hosts of even the strongest AGN can have significant star for¬ 
mation activity, unlike the trends found bv lPage et al.l d2012i) for 
radio-quiet AGN. In general, the total IR emission from the 3CR 
AGN is predominantly AGN powered, despite the frequently ac¬ 
companying strong star formation activity. 

We estimated the mass of the FIR emitting dust component, 
M dust , using 


47du S t — 


> 250/im 


D, 


K 25 Q/ 1 H! B v (250pm , / (i USt) 


( 2 ) 


where S is the flux at 250 pm rest-frame found from the 
best-fit, Dl is the luminosity distance, K 250 //m is the dust absorp¬ 
tion coefficient at 250 pm (/t 25 o^m = 4 cm 2 g -1 from the models 
of lDrainell2003l) . and B v (250pm, 7\i usl ) is the value of the Planck 
function at the corresponding rest-frame wavelength and tem¬ 
perature. Results are shown in Table [4] Given that the RGs and 
QSRs in our sample cover roughly the same redshift range and 
show similar star formation properties, it is no surprise that their 
cold dust masses are comparable as well. More interestingly, the 
masses of the cold dust component in the hosts of radio-loud 
AGN detected in at least three Herschel bands are comparable to 
those obtained for SMGs_at redshifts similar to those of our sam¬ 
ple (e.g. lSantini et al.l2010i) . The dust masses of radio-loud AGN 
provide clues to the triggering of the starburst event (and also 
that of the black hole activity, Tadhunter et al.l20~14l) . Given that 
the high-! SMGs are likely undergoing strong merger-induced 
starburst events (e.g. IKartaltepe et al.l 120121) . their similar dust 
content suggests that that FIR-luminous radio-loud AGN also 
build up their stellar mass in major gas-rich mergers. 


4.3. Median SEDs 

The best-fit SEDs for the objects detected in at least three Her¬ 
schel bands are shown in Appendix [C] The SEDs show a con¬ 
siderable range of shapes and absolute scaling, with all QSR 
(and a few RG) SEDs peaking at wavelengths around 20 pm, and 
most RG SEDs peaking at longer wavelengths. Furthermore, the 
SEDs show that the AGN-powered and star-formation-powered 
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T dus , (K) log (M dust ) (M 0 ) 



Fig. 3. Distributions of individual physical parameters obtained from the SED fits, for radio galaxies (red) and quasars (blue). The values plotted 
are those for the objects detected in at least three Herschel bands. Upper left panel: temperature of the cold dust component emitting in the FIR, 
Tdust- Upper right panel: mass of the FIR emitting cold dust component, Md USt . Lower left panel: AGN-powered IR luminosity, Lagn- Lower right 
panel: star-formation-powered IR luminosity, L S f- The vertical lines in the lower panels correspond to the average values of the FIR-detected (solid 
lines) and non-detected (dashed lines) stacked subsamples discussed in Sect. 14.41 


dust emission switch dominance typically at 50 /mi, (but this can 
happen at all wavelengths between 35 and 65 fim). Given that all 
SEDs were computed on the same rest-frame wavelength grid, 
we created median SEDs for the two types of objects to prevent 
the most extreme objects dominating the average SEDs. When 
creating the median SEDs, we refrained from applying any nor¬ 
malization, in order to preserve the absolute luminosities of the 
individual SEDs. The mediarQSEDs are shown in Fig. [5] 

The median SEDs of RGs and QSRs differ strongly at rest- 
frame 3-10 fin l, with the QSRs being a few times more luminous 
than the RGs. Such an observational difference in this wave¬ 
length regime has already been explained by lHaas et all (2 008' ) 
in the context of unification by orientation. In this scenario, the 
observed luminosity differences result from viewing the QSRs 
and RGs along different angles, such that the hot inner regions 
of the dusty torus are observed directly in the case of QSRs 
but are obscured in the case of RGs. As found by ILeinski et al.l 
(l2010i) . similar observational differences correlate with orienta¬ 
tion indicators, such as the radio core dominance. The median 
SED of RGs ca n therefore be viewed as the m edian SED of 
reddened QSRs (lHaas et al.ll2008t ILeinski et al JI201 (H) . At rest- 


4 The QSRs median SED is given only for rest-frame wavelengths 
longer than 2 fim because the emission at wavelengths shorter than 2 fim 
is dominated by the hot accretion disk, which was not included in our 
SED fitting (as explained in Appendix lAll. 


frame wavelengths between 10 and 40 fim, the median RG and 
QSR SEDs still show a considerable anisotropy, with the QSRs 
being a factor of two more luminous compared to the RGs at 
20 fim. At rest-frame FIR wavelengths (> 40 fim), however, the 
median RG and QSR SEDs appear to be remarkably similar both 
in shape and absolute scale, arguing for, on average, similar star 
formation properties for the hosts of both types of AGN. 

Figure 0 also shows the average SED for the subsample o f 
the z ~ 2 star-forming galaxies from [Kirkpatrick et al.l (120121) . 
This average SED is composed of 30 (U)LIRG galaxies, initially 
selected based on their 24 fim flux density (F^m S 100 //Jy), for 
which good multi-wavelength coverage is available throughout 
much of the IR regime. Representing the star-formation-heated 
cold d ust emission with a modified blackbody, Kirkpatric k et al.l 
(12012i ) estimate the average SFR and cold dust temperature of 
these galaxies to be 344+122 M G yr 1 and 28+2 K, respectively. 
Comparing these (U)LIRG numbers to those obtained for the 
FIR-detected 3CRs in our work, we find comparable SFRs but on 
average higher cold dust temperatures. Furthermore, Fig. [^con¬ 
firms the marked difference between the (U)LIRGs and 3CRs in 
the NIR/MIR luminosity: while the z ~ 2 star-forming galaxies 
are characterized by pronounced polycyclic aromatic hydrocar¬ 
bon (PAH) features, the powerful emission from the warm dusty 
torus completely outshines these features in the MIR SED of 
the 3CR host galaxies. Finally, the average stellar mass of the 
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Fig. 4. IR emission from star-formation-heated cold dust, Lsf, versus 
IR emission from AGN-powered dust, L AG n, for radio galaxies (filled 
red circles) and quasars (filled blue squares). Upper limits have been 
estimated as explained in the text. The large empty symbols correspond 
to the subsamples of FIR-detected RGs (circle) and QSRs (square), and 
non-detected RGs (triangle) and QSRs (diamond) discussed in Sect. 14.41 
The dashed line marks Lsf = L A gn- 



- - Median RG (3CR) 
Median QSR (3CR) 

- z~2 SF galaxies 




Rest wavelength [/mi] 


Fig. 5. Median spectral energy distributions (SEDs) for the radio galax¬ 
ies (dashed red) and quasars (dotted blue) detected in at least three Her- 
schel bands, and for z ~ 2 star-forming galaxies (solid yellow) from 
I Kirkpatrick et ah . (2012 ). Shades areas, red for radio galaxies and blue 
for quasars, correspond to the associated 16th-84th percentile ranges. 


z ~ 2 star-forming galaxies is a factor of ~ 3 smaller than that 
of the RG hosts as seen from the offset in the respective SEDs at 
around 1.6 pm. While AGN contamination to the 1.6 pm is not 
accounted for in this discussion, the very steep median SED of 
radio galaxies between 4 and 6 pm suggests that any contribution 
from hot dust to the 1.6 pm flux is negligible. AGN contamina¬ 
tion at 1.6 pm likely is present in individual cases, but does not 
appear to be crucial when discussing average properties. Over- 
all, the comparison between the 3CR and the Kirkp atrick et alJ 
d2012l) sample shows that massive galaxies in the high-z Uni¬ 
verse may be actively forming stars regardless of whether their 
supermassive black holes are accreting or not. 


Table 5. Results from the SED fitting of the stacked subsamples. 


Subsample 

log ( F A gn (Lq)) 

log (L sf (L 0 )) 

FIR-detected RGs 

12.6 

12.2 

non-detected RGs 

<12.2 

<11.7 

FIR-detected QSRs 

13.2 

12.1 

non-detected QSRs 

12.7 

<11.6 


4.4. Stacking of non-detected objects 

In addition to probing the general properties of the objects that 
are individually detected in at least three Herschel bands, we 
attempt a stacking analysis to extract an average signal for the 
objects which lack significant emission in the Herschel bands. 
For convenience, we refer to the former objects as FIR-detected, 
and to the letter as non-detected. Our aim is to discuss average 
properties for each radio-loud AGN type, therefore we split the 
non-detected objects into RGs and QSRs, respectively. In order 
to retain decent number statistics in the two subsamples, we de¬ 
cided not to further divide the RGs and QSRs in redshift bins, 
despite the considerable redshift range of our sample. In prac¬ 
tice, we take the non-detected objects to be those that are de¬ 
tected in at most one, namely the PACS 70 pm band, which given 
the redshifts of our sources often probes the peak emission from 
the dusty torus. Several RGs and QSRs were not included in the 
stacking analysis owing to the presence of potentially confusing 
sources in the corresponding maps, close to the known radio po¬ 
sition of the object in question. Furthermore, two RGs (3C 252 
and 3C 267) with a strong detection in only the PACS 70 pm 
band were not included in the stacking because they are incon¬ 
sistent with forming a single population with the objects which 
have no Herschel detections. As a result, stacking was performed 
on subsamples of 12 RGs and 6 QSRs with median redshifts 
of Zmed- 1-44 and z„,ed= 1-27, respectively. The objects entering 
the subsamples are flagged in Table [4] For comparison, we also 
selected two subsamples from the FIR-detected objects: a sub¬ 
sample with 13 RGs (Zmed- 1-34) and another one with 8 QSRs 
(Zmed= 1-52; 3C 298 and 3C 318, sources with strongly emitting 
nearby objects, were not taken into consideration). 

We stacked equal areas extracted from the individual Her¬ 
schel maps, centred on the known radio position. Photometry 
on the stacked map was performed following the same proce¬ 
dures as adopted in the case of the individually detected objects 
(Sect.H. We examined the diversity within the given subsample 
by bootstrapping with 1000 realizations. In practice, from the 
original subsamples identified above, for each bootstrapping re¬ 
alization we selected a random subsample (with the same num¬ 
ber of objects, allowing for repetitions), stacked the Herschel 
maps, and performed the photometry. The centroid and the dis¬ 
persion of the resulting distribution were taken to be the mean 
flux density of the on-source stack and its associated uncertainty. 
Additionally, we stacked random positions in order to inspect the 
overall significance of the on-source stacked signal. If the mean 
value of the on-source stack distribution was at least three times 
larger than the mean value of the corresponding background 
stack distribution, we concluded that the on-source signal is sig¬ 
nificant. In cases of a non-significant signal, we took three times 
the mean value of the background stack distribution to be our 
on-source stack upper limit value. 

The mean stacks for the four different subsamples selected 
in our study are presented in Fig. [6] The non-detected quasars 
have a significant stacked signal in only the PACS 70 pm band. 
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Fig. 6. Stacked images in all Herschel bands, for the four subsamples discussed here. From left to right: PACS 70 /mi, PACS 160 /mi, SPIRE 
250 /mi, SPIRE 350 /on, and SPIRE 500 /mi bands, respectively. From top to bottom: FIR-detected radio galaxies, non-detected radio galaxies, 
FIR-detected quasars, and non-detected quasars. Each stacked image shown here has dimensions of 2x2 arcmin and is centred on the known radio 
positions of the objects entering the stack. 



Fig. 7. Spectral energy distributions of stacked subsamples selected 
from the hosts of the 3CR sources studies here. The four subsamples 
presented are: FIR-detected radio galaxies (red circles), non-detected 
radio galaxies (red triangles), FIR-detected quasars (blue squares), and 
non-detected quasars (blue diamonds). Arrows mark 3cr upper limits as 
explained in the text. 


whereas the non-detected radio galaxies have no significant sig¬ 
nal in any of the Herschel bands. On the other hand, the FIR- 


detected quasars and radio galaxies have significant stacked sig¬ 
nals in all, and in all but the SPIRE 500 //m band, respectively. 
In order to study the full IR SED of a stacked subsample, we 
calculated mean flux densities in the Spitzer bands, where vir¬ 
tually all objects have been individually detected. All obtained 
average fluxes are in the observed frame: these were transferred 
to the rest-frame, using the median redshift of the corresponding 
subsample. The SEDs of the stacked subsamples are shown in 
Fig. ID 

The 1-10 pm SEDs of radio galaxies are composed of emis¬ 
sion from an evolved stellar population and from a heavily red¬ 
dened AGN dHaas et al.l 12008 ). At the rest-frame wavelengths 
(given the comparable median redshifts of the various sub¬ 
samples) probed by the two shortest wavelength Spitzer bands 
(< 2 pm), the FIR-detected and non-detected subsamples of 
RGs are similar. Thus, it is likely that the stellar masses of 
the host galaxies for the two RG subsamples are, on average, 
similar as well. Identifying differences in the host galaxy stel¬ 
lar masses requires additional NIR photometry, and is outside 
the scope of this work. At wavelengths between 2 and 10 pm, 
the SEDs of the two RG subsamples differ from each other. At 
these wavelengths, AGN-heated hot dust emission appears to be 
present/absent in the SED of the FIR-detected/non-detected sub¬ 
sample, respectively, reflecting different levels of dust obscura¬ 
tion, as inferred from the diversity in the individual MIR SEDs 
of RGs dHaas et al J 12008 ). Particularly interesting is the non- 
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significant stacked signal for non-detected RGs at rest-frame 
30 pm, which could indicate that the two subsamples have dif¬ 
ferent intrinsic AGN luminosities. Indeed, fitting the SEDs fol¬ 
lowing the routine taken in Sect. [3] gives (at least) a factor of 
two difference between the estimated Lagn in the two RGs sub¬ 
samples (Table |5] Fig.Q. A clear difference between the SEDs 
of the two subsamples of RGs is observed at rest-frame wave¬ 
lengths longer than 60 pm, suggesting marked differences in the 
star formation properties of the two subsamples. This results in 
(at least) a factor of three difference between the estimated SFRs 
in the two RGs subsamples (Table [3 Fig. []}• The SEDs of the 
two QSR subsamples generally follow the same trends as those 
observed for the RG subsamples. While the non-detected QSRs 
could overall be at the fainter end of the QSR population, the 
most striking difference in the QSR SEDs is at FIR wavelengths, 
with the FIR-detected QSRs having SFRs (at least) a factor of 
three higher than the non-detected QSRs. While more analysis 
is required to pin down the potential differences among the sub¬ 
samples in wavelengths other than FIR, it is beyond doubt that 
the formation of new stars can be prodigious in some, but modest 
or weak in other radio-loud AGN hosts. 


5. Discussion 

5.1. Star formation in hosts of powerful AGN 

Starbursts powering the FIR emission in some hosts of radio- 
loud AGN_have been found at low-to-intermediate redshifts (e.g. 
iDicken et al.l l2Q10h . In the high-’ Universe, high levels of star 
formation in some radio-loud AGN have been estimated using 
several observational indicators. These include the usage of rest- 
frame UV spectroscopy to detect huge Fya halos surrounding 
the AGN (e.g. IVillar-Martm et alii 1999!) . submm photometry to 
probe the cool dust and molecular gas conten t of AGN hosts (e.g. 
lArchibald et al.l 1200 It iReuland et al.l 120041) . an d Spitzer MIR 
spectr oscopy to detect strong PAH features (e.g. lRawlings et all 
120131 ). The latter, however, are seen only in rare cases because 
the AGN-powered hot/warm dust emission in powerful radio- 
loud AGN hosts usually outshines the PAH emission. 

UV/visible data have been used to infer SFRs for several 
3CR radio galaxies studied in this work (Chambers & Chariot) 
1990). These include 3C 065,3C 068.2,3C 266,3C 267,3C 324, 
3C 356 and 3C 368. The ratio between the SFR obtained us¬ 
ing the approach taken in this work and the UV/visible SFR 
ranges between 4 and 40, therefore we conclude that the star 
formation is often strongly obscured in the UV/visible. This 
comparison demonstrates that rest-frame FIR data are crucial 
in quantifying the dust-enshrouded star formation in the hosts 
of these sources, especially the ones hosting quasars because 
the emission from the accretion disk complicates the SFRs esti¬ 
mates from UV/visible data. Using Herschel data. ISevmour et akl 
d201 ill found mean SFRs for 1.2 < z < 3.0 radio-selected AGN 
to range between 80 and 581 M 0 yr 1 . Similarly. iDrouart et al.l 
(120141) studied a sample of ~ 70 1 < z < 5 radio galaxies, es¬ 
timating SFRs of a few hundred to a few thousand solar masses 
per year. While the authors suggest that the highest values re¬ 
ported in their study are likely overestimated and should be 
treated as upper limits, the overall idea that the hosts of high- 
Z radio-loud AGN can be prodigious star-formers is consistent 
with our findings for the hosts of the powerful 3CR objects. A 
crucial point in these analyses is that the high star formation lu¬ 
minosities obtained are coeval with the growth of the SMBHs 
residing in the nuclei of the host galaxies. As such, these re¬ 
sults argue for a scenario whereby on average the SMBH has not 


quenched the star formation in the host galaxy, which is at odds 
with results presented by Page et akl d2012l) . A series of coeval 
episodes of strong star formation and black hole activity may 
have formed the massive host galaxies and their massive black 
holes (see B12). 

The SFRs estimated for the FIR-detected radio-loud AGN 
are comparable to those obtained for SMGs at similar redshifts 
(e.g. iMagnelli et akil2012l) . Other parameters, such as the tem¬ 
peratures and masses of the cold dust component and total stel- 
lar masses agree w it h each other, at least on average, as well 
(ISantini et all l20fot iMichalowski et al] 120121 : ISwinbank et al.l 
2014). It is widely thought that high-z SMGs form stars in star- 
burst events, induced as a result of a variety of different processes 
including mergers and tidal interactions. Radio-loud AGN are 
known to be at the centres of over-densities in the high -z Uni¬ 
verse (e.g. IVenemans et akll2007l : IWvlezalek et alJl2013h . Thus, 
the merger scenario appears to be an attractive way of produc¬ 
ing at least some of the extremely high SFRs estimated in our 
study. A similar conclusion has been drawn in a study of a few 
high-z_radio galaxies using both Herschel and CO observations 


__ 

(flvison et al.li2012t ). 

Jet-induced star formation, also known as positive feed back, 
is yet another possible way to trigger high SFRs (e.g. lDev et al.1 
119971 : Ivan Breugel et alJfl998h in hosts of radio-loud AGN. In 
this scenario, the outgoing radio jet shocks the surrounding in¬ 
terstellar material, which subsequently cools down to form new 
stars. Recently, based on UV-to-submm templates built with the 
evolutionary code PEGASE.3. iRocca-Volmerange et all (120131) 
showed that the star formation timescales from stellar popula¬ 
tion synthesis agree well with the ages of the radio episodes for 
two high-z radio galaxies. Similar studies of larger samples are 
required to better understand the details of jet-induced star for¬ 
mation. 

Despite evidence for strong star formation in the 3CR hosts 
detected in at least three Herschel bands, ~ 60% of 3CR sources 
are detected in fewer than three Herschel bands, arguing for 
significantly lower star formation activity in these hosts. While 
the upper limits in the SPIRE bands still allow energetically 
significant SFRs of up to 300 M e yr~ 1 in some of these hosts, 
the overwhelming majority have SFRs of at most 100 M 0 yr 1 . 
This raises an important question: what leads to the significantly 
lower star formation activity in these hosts? If the star forma¬ 
tion activity is indeed merger-induced, one possible answer to 
the above question is that occurrences of minor- and/or gas-poor 
mergers are likely to lead to relatively low levels of star forma¬ 
tion activity. The relevance of the merger scenario in the con¬ 
text of the triggering of starbursts, but also o f the AGN activity, 
has been discussed by Tadhunter et alJ (1201 ll) . They find that the 
strongly starbursting radio galaxies in their intermediate-redshift 
sample have optical morphological features consistent with the 
idea that they are triggered in major mergers. Another possi¬ 
bility is that we are observing some hosts after the star forma¬ 
tion has been quenched (e.g. lFarrah et al.ll2012h . If the AGN ac¬ 
tivity is responsible for the quenching, then we argue that this 
negative feedback is not universal, even if it acts over a very 
short timescale. First, given the strong star formation in the 3CR 
hosts detected in at least three Herschel bands, it is unlikely that 
the quenching of star formation occurs before the radio-loud 
AGN phase. Second, quenching taking place during the radio- 
loud phase is not supported with the finding in Sect. l4.4l that the 
< 2 pm rest-frame stacked SEDs of the two RG subsamples are 
similar. In particular, if star formation in the non-detected RGs 
is quenched, then they are expected to be brighter than the FIR- 
detected RGs because the star formation in the latter is heav- 
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Fig. 8. Comparison of star formation rates of AGN and non-AGN host 
galaxies. Plotted are the objects detected in at least three Herschel bands 
(dashed green). A comparison sample of objects (solid magenta) is se¬ 
lected from some of the deepes t_Herschel data, from PACS-GOODS 
South field (Rodiehiero et al. 2011). The vertical lines correspond to 
the average values of the FIR-detected (solid lines) and non-detected 
(dashed lines) stacked subsamples, discussed in Sect. 14.41 


ily dust obscured. Third, because the radio-loud phase probably 
marks the end of the AGN phase, it is unlikely that the quenching 
occurs after the radio-loud phase. To conclude, while the hosts of 
the quasar-mode radio-loud AGN studied here may have a wide 
range of star formation rates, it is unlikely that the triggering or 
quenching of their star formation activity is associated with a 
uniform scenario met in all objects. 

5.2. RG and QSR unification in the FIR 

Unification the ories of radio-loud AGN (Barthell Il989t 
lAnton ucci 1993) ascribe observed differences in the properties 
of radio galaxies and quasars solely to orientation effects. Us¬ 
ing Spitzer photometric observations of the sample of 3CR ob¬ 
jects studied here, (iHaas et all2008l) (see also the median SEDs 
from our current work) found the mean 1-10 pm rest-frame ra¬ 
dio galaxy SED to be consistent with a sum of an underly¬ 
ing host and a heavily obscured quasar. Unification among the 
high-z 3CR objects was further corroborated with Spitzer spec- 
troscopic dLeipski et al.l l2010) and Chandra X-ray observations 
(IWilkes et al.l2013l) . Our results reveal that the optically thin (i.e. 
isotropic) FIR emission is similar for radio galaxies and quasars, 
thus in line with the predictions of radio-loud AGN unification 
by orientation. 


5.3. AGN versus non-AGN host galaxies 


It is important to check how the SFRs obtained for the FIR- 
detected objects in this work compare to those for the non- 
AGN galaxy population. SFRs in star-forming galaxies are 
strongly correlated with both stellar mass and redshift (e.g. 
lElbaz et al.ll2007l) . therefore, the comparison must be made with 
a sample having a comparable range of stellar mass and red- 
shift. About a dozen 3CR radio galaxies studied in this work 
have stellar mass estimates primarily based on SED-fitting of 
visibl e/NIR dBest et al. [l998) or MIR broad-band photome¬ 
try (ISevmour et al. 2007tlDe Breuck et al.ll2010h . The estimated 
stellar masses range between 1.5 x 10 n M 0 and 6 x 10 11 M 0 , and 


are not e xpected to be a fu nction of the redshift/luminosity of the 
sources dDe Breuck et al .ioToh. Measurements of stellar masses 
in high-’ quasar hosts are problematic because the strong con¬ 
tinuum emission from the accretion disk often outshines their 
host galaxies. However, assuming that unification of powerful 
radio galaxies and quasars holds, the masses of quasars and ra¬ 
dio galaxies hosts are often taken to be (at least on average) sim¬ 
ilar (e.g. lMcLure et alj2006i) . Thus, for the discussion below, we 
assume that the stellar masses of our 3CR hosts range between 
1.5 x 10 n M o and 6 x 10 n M o 

As a control sample of non-AGN galaxies, we selected 
1.5 < z < 2.5 galaxies within the stellar mass range in¬ 
dicated above, whose star-forming properties were estimated 
from deep PACS data of the GOODS South field (cyan points 
in Fig.l from iRodighiero et al.1 l20l7l) . The stellar masses of 
the se galaxies were e stima ted from SED-fitting as explained 
by iRodighiero et al. (2 01(1) . While the majority of the se¬ 
lected galaxies lie on the main sequence of star-forming galax¬ 
ies, a fraction are located above it, most likely character¬ 
ized by the starbursting nature of the ongoing star formation 
(iRodighiero et alJ 1201 ill . Figure [8] shows the SFR histograms 
of the 3CRs and the selected control sample. It is clear that 
the 3CR FIR-detected objects have, on average, SFRs compa¬ 
rable to those of their equally massive non-AGN counterparts. 
The majority of FIR-detected objects are thus also located near 
the main sequence of star-forming galaxies, similarly to what 
has been found in deep Herschelf PACS studies of less powerful, 
X-rav selected, high -z AGN (iMullanev et alJl2012b ISantini et al.l 
1201 2t [Rosario et afl 1 201 3l) . In comparison to the FIR-detected 
objects, the non-detected objects have similar stellar mass but 
significantly lower SFRs, placing them below the main sequence 
of star-forming galaxies. We cannot exclude the possibility that 
the star formation activity in some of these objects has been 
quenched. Better estimates of the stellar masses of our high-z 
3CR sources will allow a more robust statistical study of their 
exact location with respect to the main sequence of star-forming 
galaxies. 

5.4. Subgalactic versus supergalactic radio sources 

The radio morphologies of radio-loud AGN present a unique op¬ 
portunity to estimate the duration of the AGN episode by as¬ 
suming a typical speed of radio jet expansion of 10%-20% of 
the speed of light. The projected radio sizes of our high-z 3CR 
sources, measured lobe-to-lobe, have been measured from high- 
resolution radio images. Based on their projected radio sizes, we 
divided the 3CR sources into two groups, subgalactic (< 30 kpc) 
and supergalactic radio (> 30 kpc) sources. The subgalactic 
sources, typically contained within their host galaxies, account 
for 25% of our high-z 3CR sample. Figure [9] shows the objects’ 
estimated SFRs as a function of their projected radio size. Both 
the subgalactic and supergalactic FIR-detected 3CR hosts have 
comparable SFRs. In contrast, the majority of the non-detected 
objects have larger, i.e. older radio sources. As such, the ratio of 
FIR-detected versus non-detected objects appears to be a func¬ 
tion of projected radio size, changing from 1.3 for subgalactic to 
0.5 for supergalactic sources. 

There are at least two different effects contributing to the 
findings presented in Fig. [9] First, the finding that many qui¬ 
escent galaxies turn up when the radio sources are large is con¬ 
sistent with the fact that star formation depends heavily on the 
availability of cold gas. Indeed, the process of exhausting the 
available fuel for star formation has timescales similar to the age 
of large radio sources. Second, our result that many small radio 
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Fig. 9. Estimated star formation rate as a function of projected radio size 
(measured lobe to lobe) for the 3CR galaxies (red circles) and quasars 
(blue squares). The dashed line indicates the value taken to differenti¬ 
ate between subgalactic (< 30 kpc) and supergalactic (> 30 kpc) radio 
sources, respectively. Arrows indicate star formation rate upper limits 
for objects with fewer than three Herschel detections. 


sources are hosted by strongly star-forming galaxies is consistent 


with the observational/theoretical finding that radi 
duce bursts of star formation (positive feedback, e 

j jets may in- 

g. Silk 

2013) 

in the hosts of high-z (Dev et al. 1997; van Breugel et al. 

1998) 

and low-z radio-loud AGN (Tadhunteret al. 2011; 

Dicken et al. 


l2012l) . lBest et alJ (l996) found that smaller radio sources show 
stronger alignment effect providing further support for the in¬ 
cidence of positive feedback in hosts of small radio sources. A 
step forward in probing the incidence of positive feedback within 
the hosts of 3CR sources may be achieved by correlating the 
ages of the young stellar components and those of the current ra¬ 
dio episode, similarly to the study by iRocca-Volmerange et al] 
( 2013 (1). An issue which complicates the overall picture, how¬ 
ever, is whether the link between the smaller radio sources and 
star formation activity is a consequence of an observational bias. 
Namely, as pointed out by iTadhunter et alJ (1201 ill , the interac¬ 
tion between the jet and the host ISM in subgalactic sources may 
boost the radio emission, leading to a preferential selection of 
such sources in flux-limited samples, like the 3CR studied here. 

5.5. Model fit limitations 

Our estimates of the cold dust temperature may suffer from a 
bias due to the known degeneracy between the cold dust temper¬ 
ature and the dust emissivity index, J3: a lower fixed fi value will 
lead to a higher dust temperature. Constraining the emissivity in¬ 
dex within our sample is a difficult task because (1) the peak of 
the dust emission is not well isolated in the SED, and (2) only a 
few data points probe the Rayleigh-Jeans tail of the dust emis¬ 
sion. Because of this, and to limit the number of free parameters 
throughout the fitting, we fixed the emissivity index to a value of 

1.6. However, the detections in all Herschel bands, in addition to 
the availability of submm data, allow /3 to be constrained for a 
few selected objects. The ft values for these five objects (flagged 
in Table Q range from 1 to 2.3, with increased/decreased esti- 

5 The alignment effect is the co-spatial extent of radio and 
UV/visible/NIR emission in radio-loud AGN, partly due to the inter- 
action between the jet and the interstellar matter of the host galaxy 
dMcCarthvll 1 9931: iMilev & D e Breuck 2008). 


mated cold dust temperatures for lower/higher p values, respec¬ 
tively. Nevertheless, the star formation luminosities, and conse¬ 
quently SFRs, remain the same (within 10%). 

While we attribute the FIR emission to star formation on the 
scale of the host galaxy, we stress that the FIR emission consid¬ 
ered is the excess emissio n after the AGN-powered to rus emis¬ 
sion (using the models of Honig & Kishimoto 2010) has been 
accounted for. From our SEDs, we concluded that star forma¬ 
tion is the dominant process yielding emission at_ rest-frame 
wavelengths longer than A > 50 pm (see also iFeipski et al.l 
120131) . Nevertheless, should future AGN models demonstrate 
that the AGN-powered emission continues to dominate at rest- 
frame wavelengths longer than A ~ 50 pm, then our SFR esti¬ 
mates are likely to be upper limits. 

6. Conclusions 

We present Herschel photometry of the complete sample of z > 1 
3CR radio galaxies and quasars. The 3CR sample is a flux- 
limited sample, consisting of some of the most powerful radio- 
loud AGN accreting in quasar-mode. Combining the Herschel 
photometry with available Spitzer data, we performed a full IR 
SED analysis, separating the contribution from the AGN and 
from the star formation activity in the host galaxy. We summa¬ 
rize our findings below: 

1. About 40% of the studied objects have robust PACS 
and SPIRE detections, translating into UFIRG-like star 
formation luminosities, i.e. of the order of SFRs of hundreds 
of solar masses per year. Such prodigious levels of star 
formation have recently also been inferred for other high-z 
radio galaxies (e.g. IDrouart et al.l 12014)) . Merger induced 
and/or jet triggered star formation activity are both possible 
mechanisms leading to the SFRs obtained for these objects. 

2. The SFRs of the FIR-detected objects are comparable to 
those of mass-matched, non-AGN galaxies, selected from 
deep Herschel surveys. There is no clear evidence that 
the star formation has been quenched in the hosts of the 
FIR-detected objects. 

3. The total IR (1-1000 pm) emission from the high-z 3CR 
radio sources is predominantly powered by the AGN, despite 
the frequently strong starburst activity coeval with the AGN 
episode. Furthermore, no strong correlation between the 
AGN- and star formation powered IR luminosities is found. 

4. The median SEDs of the FIR-detected objects show that 
RGs and QSRs are quite different in the MIR, but remark¬ 
ably similar in the FIR. Thus, while the MIR emission is 
anisotropic, the FIR emission is isotropic and optically thin. 
These findings are consistent with the orientation-based 
unification of radio-loud AGN. 

5. Splitting the sample into subgalactic (< 30 kpc) and 
supergalactic (> 30 kpc) radio sources, the fraction of 
Herschel detected objects is a function of the projected radio 
size of the sources. In particular, the hosts of subgalactic 
radio sources are more likely to be detected by Herschel, 
arguing for a possible link between radio size and star 
formation activity, i.e. jet-induced star formation (positive 
feedback), or for fading of star formation in mature AGN. 
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6. Stacking of the Herschel non-detected objects reveals a 
class of MIR/FIR faint objects. While ongoing star formation 
episodes (at significantly lower levels than those discussed 
above) cannot be ruled out, star formation has largely ceased 
in the hosts of these objects. As such, the radio-selected, 
high-;; 3CR hosts appear to be a heterogeneous mixture of 
MIR/FIR bright and faint objects. 

Upcoming instruments with better sensitivity/resolution, 
such as ALMA, will likely help us pinpoint the exact location 
of the ongoing star formation in high-z galaxies. This will lead 
to a further understanding of the interplay between the AGN and 
star formation activity within high-z AGN hosts. 
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Table 1. Summary of selected properties of the high-z 3CR sample studied in this work. 


Name 

(1) 

Type 

(2) 

z 

(3) 

RA (J2000) 
(4) 

Dec (J2000) 
(5) 

log (L 178MH z (W Hz- 1 )) 
(6) 

PACS OBSIDs 

(7) 

SPIRE OBSIDs 

(8) 

3C 002 

QSR 

1.04 

00:06:22.58 

-00:04:24.69 

29.0 

1342221796/1342221797 

1342212374 

3C 009 

QSR 

2.01 

00:20:25.21 

+ 15:40:54.59 

29.8 

1342222444/1342222445 

1342213198 

3C013 

RG 

1.35 

00:34:14.55 

+39:24:16.65 

29.2 

1342223179/1342223180 

1342213491 

3C014 

QSR 

1.47 

00:36:06.44 

+ 18:37:59.23 

29.2 

1342222429/1342222428 

1342213196 

3C 036 

RG 

1.30 

01:17:59.48 

+45:36:21.75 

29.0 

1342223508/1342223509 

1342203613 

3C043 

QSR 

1.47 

01:29:59.80 

+23:38:20.28 

29.2 

1342223506/1342223507 

1342213488 

3C 065° 

RG 

1.18 

02:23:43.19 

+40:00:52.45 

29.1 

1342238005/1342238006 

1342239821 

3C 068.1 

QSR 

1.24 

02:32:28.87 

+34:23:46.79 

29.1 

1342223870/1342223871 

1342226628 

3C 068.2 

RG 

1.58 

02:34:23.85 

+31:34:17.46 

29.2 

1342223866/1342223867 

1342224971 

3C 119 

RG 

1.02 

04:32:36.50 

+41:38:28.44 

29.0 

1342227975/1342227976 

1342216924 

3C 124 

RG 

1.08 

04:41:59.10 

+01:21:01.91 

28.9 

1342226718/1342226719 

1342216939 

3C 173 

RG 

1.03 

07:02:20.58 

+37:57:23.50 

28.7 

1342219418/1342219419 

1342206177 

3C 181 

QSR 

1.38 

07:28:10.30 

+ 14:37:36.24 

29.3 

1342220573/1342220574 

1342204852 

3C 186 

QSR 

1.06 

07:44:17.45 

+37:53:17.15 

29.0 

1342220127/1342220128 

1342206178 

3C 190 

QSR 

1.20 

08:01:33.55 

+ 14:14:42.94 

29.1 

1342205262/1342205263 

1342205052 

3C 191 

QSR 

1.96 

08:04:47.97 

+ 10:15:23.69 

29.6 

1342220655/1342220656 

1342205072 

3C 194 

RG 

1.18 

08:10:03.61 

+42:28:04.31 

28.9 

1342220123/1342220124 

1342206180 

3C 204 

QSR 

1.11 

08:37:44.95 

+65:13:34.92 

28.9 

1342220115/1342220116 

1342206190 

3C 205 

QSR 

1.53 

08:39:06.45 

+57:54:17.12 

29.3 

1342220117/1342220118 

1342206188 

3C 208.0 

QSR 

1.11 

08:53:08.60 

+ 13:52:54.98 

29.1 

1342220790/1342220791 

1342206221 

3C 208.1 

RG 

1.02 

08:54:39.28 

+ 14:05:52.56 

28.7 

1342220788/1342220789 

1342206220 

3C210 

RG 

1.17 

08:58:09.96 

+27:50:51.57 

28.9 

1342220797/1342220796 

1342230768 

3C212 

QSR 

1.05 

08:58:41.49 

+ 14:09:43.97 

29.0 

1342220786/1342220787 

1342206219 

3C 220.2 

QSR 

1.16 

09:30:33.47 

+36:01:24.17 

28.8 

1342220798/1342220799 

1342222125 

3C 222 

RG 

1.34 

09:36:32.01 

+04:22:10.30 

29.1 

1342221142/1342221143 

1342210521 

3C 225A 

RG 

1.56 

09:42:08.48 

+ 13:51:54.23 

29.1 

1342221258/1342221259 

1342210518 

3C 230 

RG 

1.49 

09:51:58.82 

-00:01:27.23 

29.5 

1342221136/1342221137 

1342210520 

3C 238 

RG 

1.40 

10:11:00.37 

+06:24:39.72 

29.3 

1342221144/1342221145 

1342210519 

3C 239" 

RG 

1.78 

10:11:45.41 

+46:28:19.75 

29.5 

1342231241/1342231242 

1342230739 

3C 241 

RG 

1.62 

10:21:54.52 

+21:59:30.71 

29.3 

1342221152/1342221153 

1342198253 

3C 245 

QSR 

1.03 

10:42:44.60 

+ 12:03:31.26 

28.9 

1342221264/1342221265 

1342210516 

3C 249 

RG 

1.55 

11:02:03.84 

-01:16:17.39 

29.5 

1342221853/1342221854 

1342198569 

3C 250 

RG 

1.26 

11:08:52.12 

+25:00:54.61 

29.1 

1342221154/1342221155 

1342210509 

3C 252 

RG 

1.10 

11:11:32.99 

+35:40:41.64 

28.9 

1342221160/1342221161 

1342210508 

3C 255 

RG 

1.36 

11:19:25.23 

-03:02:51.50 

29.2 

1342221851/1342221852 

1342210515 

3C 256 

RG 

1.82 

11:20:43.02 

+23:27:55.22 

29.4 

1342221262/1342221263 

1342210510 

3C 257 

RG 

2.47 

11:23:09.17 

+05:30:19.47 

29.7 

1342221966/1342221967 

1342210514 

3C 266" 

RG 

1.27 

11:45:43.36 

+49:46:08.24 

29.1 

1342222695/1342222696 

1342222663 

3C 267 

RG 

1.14 

11:49:56.56 

+ 12:47:19.07 

29.1 

1342222448/1342222449 

1342200236 

3C 268.4 

QSR 

1.40 

12:09:13.61 

+43:39:20.96 

29.1 

1342221162/1342221163 

1342210501 

3C 270.1 

QSR 

1.52 

12:20:33.87 

+33:43:12.05 

29.3 

1342221952/1342221953 

1342200238 

3C 280.1 

QSR 

1.66 

13:00:33.30 

+40:09:07.72 

29.3 

1342212393/1342212394 

1342210498 

3C 294" 

RG 

1.78 

14:06:53.20 

+34:11:21.10 

29.4 

1342211098/1342211099 

1342206200 

3C 297 

RG 

1.41 

14:17:23.99 

-04:00:47.54 

29.1 

1342223834/1342223835 

1342203577 

3C 298 

QSR 

1.44 

14:19:08.18 

+06:28:34.80 

29.8 

1342223664/1342223665 

1342213464 

3C 305.1 

RG 

1.13 

14:47:09.56 

+76:56:21.80 

29.0 

1342220952/1342220953 

1342206193 

3C318 

QSR 

1.57 

15:20:05.44 

+20:16:05.76 

29.3 

1342223844/1342223845 

1342204107 

3C 322 

RG 

1.68 

15:35:01.23 

+55:36:52.87 

29.3 

1342199131/1342199132 

1342206196 

3C 324 

RG 

1.21 

15:49:48.89 

+21:25:38.06 

29.2 

1342202562/1342202563 

1342213461 
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Table 1. continued. 


Name 

(1) 

Type 

(2) 

z 

(3) 

RA (J2000) 
(4) 

Dec (J2000) 
(5) 

log (L 178MH z (W Hz- 1 )) 
(6) 

PACS OBSIDs 

(7) 

SPIRE OBSIDs 

(8) 

3C 325 

QSR 

1.13 

15:49:58.42 

+62:41:21.66 

29.1 

1342219034/1342219035 

1342206195 

3C 326.1 

RG 

1.83 

15:56:10.06 

+20:04:20.44 

29.3 

1342224482/1342224483 

1342213462 

3C 356 

RG 

1.08 

17:24:19.04 

+50:57:40.14 

28.9 

1342219036/1342219037 

1342206197 

3C 368 

RG 

1.13 

18:05:06.45 

+ 11:01:35.06 

29.0 

1342216599/1342216600 

1342216954 

3C418 

QSR 

1.69 

20:38:37.03 

+51:19:12.66 

29.4 

1342219032/1342219033 

1342210542 

3C 432 

QSR 

1.80 

21:22:46.32 

+ 17:04:37.95 

29.4 

1342211499/1342211500 

1342210541 

3C 437 

RG 

1.48 

21:47:25.10 

+ 15:20:37.49 

29.3 

1342211497/1342211498 

1342210540 

3C 454.0 

QSR 

1.76 

22:51:34.73 

+ 18:48:40.12 

29.4 

1342210949/1342210950 

1342210539 

3C 454.1 

RG 

1.84 

22:50:32.93 

+71:29:19.18 

29.4 

1342211436/1342211437 

1342212365 

3C 469.1 

RG 

1.34 

23:55:23.32 

+79:55:19.60 

29.1 

1342221170/1342221171 

1342220543 

3C 470" 

RG 

1.65 

23:58:35.89 

+44:04:45.55 

29.3 

1342237858/1342237859 

1342236248 

4C 13.66 

RG 

1.45 

18:01:38.95 

+ 13:51:23.85 

29.2 

1342216597/1342216598 

1342216956 

4C 16.49 

QSR 

1.30 

17:34:42.61 

+ 16:00:31.21 

29.1 

1342216595/1342216596 

1342216955 


Notes. (1) Name of object; (2) AGN type; (3) Redshift; (4) Right Ascension; (5) Declination; (6) Log of 178 MHz luminosity (in the observer 
frame); (7) PACS ObsID; (8) SPIRE ObsID 

(a> Herschel observations taken from OTl_nseymour_l (PI: Seymour) 


Article number, page 15 of[30] 






Article number, page 16 of®] 


Table 2 . Herschel and Spitzer photometry of the 3CR objects studied in this work. Photometric uncertainties are lcr values, and upper limits are 3<x values. 


Object 

( 1 ) 

OJy) 

( 2 ) 

F 4.5//m 

(>Jy) 

(3) 

F 5 . 8 // 1 T 1 

(>Jy) 

(4) 

Fg.O/m 

OJy) 

(5) 

Fl 6 /im 

OJy) 

( 6 ) 

F24/jm 

GuJy) 

(7) 

F7(Vm 

(mJy) 

( 8 ) 

F 160/im 
(mJy) 

(9) 

F 250 /J 111 

(mJy) 

( 10 ) 

F 35 O /1111 

(mJy) 

(11) 

F500/jm 

(mJy) 

(12) 

3C 002 

283+42 

330+50 

530+80 

809+121 

1550+233 

2970+446 

10.0+2.6 

21.9+4.4 

16.1+5.1 

19.6+4.8 

< 17.0 

3C 009 

884+133 

1080+162 

1590+239 

2220+333 

3330+500 

3470+520 

13.7+2.6 

18.1+4.9 

< 15.4 

< 13.8 

< 17.2 

3C013 

133+20 

133+20 

147+22 

283+42 

375+56 

2060+309 

24.6+2.1 

30.7+5.8 

< 18.8 

< 15.3 

<20.7 

3C014 

1040+156 

1710+257 

2740+411 

4150+623 

7070+1061 

10300+1545 

20.9+2.6 

21.9+6.8 

20.0+6.3 

< 15.4 

< 19.8 

3C 036 

163+24 

205+31 

256+38 

360+54 

560+84 

874+131 

<4.2 

< 10.6 

< 11.2 

< 11.8 

< 14.2 

3C043 

193+29 

270+41 

356+53 

445+67 

1010+152 

1610+242 

b 

b 

< 14.1 

< 17.9 

< 15.8 

3C 065 

202+30 

233+35 

299+45 

418+63 

798+120 

1700+255 

< 6.7 

< 12.1 

<20.8 

< 18.1 

<20.2 

3C 068.1 

967+145 

1430+215 

2040+306 

2780+417 

3800+570 

7760+1164 

22.7+2.3 

< 17.8 

< 18.2 

< 15.9 

<20.1 

3C 068.2 

105+16 

129+19 

137+21 

112+17 

1340+201 

1170+176 

27.4+2.6 

39.6+5.9 

42.0+7.2 

38.7+7.0 

29.0+7.2 

3C 119 

802+120 

878+132 

1280+192 

1850+278 

4820+723 

8260+1239 

24.9+2.2 

32.6+8.9 

28.9+15.4“ 

<56.1 

<44.1 

3C 124 

144+22 

120+18 

188+28 

310+47 

1840+276 

3560+534 

34.2+2.1 

55.7+6.7 

52.1+7.3 

31.5+7.9 

<24.3 

3C 173 

163+24 

172+26 

197+30 

227+34 

374+56 

710+107 

6.9+1.6 

< 10.2 

< 11.6 

< 13.0 

< 16.5 

3C 181 

348+52 

485+73 

722+108 

1110+167 

2180+327 

4260+639 

12.2+2.3 

< 13.6 

<20.4 

< 16.5 

<21.8 

3C 186 

791+119 

1020+153 

1410+212 

1960+294 

3660+549 

6660+999 

18.9+2.6 

< 18.0 

< 13.7 

< 17.1 

<21.0 

3C 190 

739+111 

908+136 

1290+194 

1740+261 

3310+497 

6690+1004 

46.1+2.5 

72.5+4.8 

74.1+6.5 

54.1+4.7 

<20.9 

3C 191 

333+50 

399+60 

655+98 

1010+152 

2270+341 

3810+572 

26.4+2.7 

21.9+4.9 

< 18.9 

< 17.1 

< 18.1 

3C 194 

201+30 

176+26 

164+25 

208+31 

509+76 

885+133 

<4.6 

<9.3 

<20.2 

< 18.9 

<20.7 

3C 204 

917+138 

1250+188 

1920+288 

2540+381 

4730+710 

7360+1104 

<6.1 

< 13.5 

< 17.9 

< 14.5 

< 18.2 

3C 205 

1460+219 

2080+312 

2920+438 

4090+614 

7320+1098 

12800+1920 

62.7+2.6 

66.2+6.8 

56.2+6.8 

31.6+4.8 

19.9+6.7“ 

3C 208.0 

660+99 

803+120 

1160+174 

1620+243 

2980+447 

5870+881 

<7.4 

< 16.7 

<20.5 

<21.7 

<20.3 

3C 208.1 

331+50 

430+65 

656+98 

954+143 

1360+204 

2110+317 

<8.1 

< 15.7 

< 17.1 

< 16.0 

<20.2 

3C210 

256+38 

336+50 

489+73 

1090+164 

3410+512 

4430+665 

31.6+2.4 

56.0+4.0 

b 

b 

b 

3C212 

925+139 

1430+215 

2340+351 

3400+510 

6710+1007 

10800+1620 

16.6+2.6 

< 16.7 

<38.5 

<56.0 

< 39.9 

3C 220.2 

592+89 

870+131 

1330+200 

2000+300 

4150+623 

6720+1008 

26.6+2.2 

22.5+5.3 

< 14.1 

< 14.6 

< 15.7 

3C 222 

83+12 

91 + 14 

73+11 

65+10 

331+50 

229+34 

14.3+1.9 

50.8+4.7 

48.4+4.8 

50.9+3.5 

28.7+5.4 

3C 225A 

47+7 

49+7 

71 + 11 

108+16 

321+48 

< 1070 

<7.8 

<21.4 

<21.3 

< 16.3 

<22.7 

3C 230 

1040+156 

672+101 

438+66 

317+48 

1150+173 

1560+234 

11.9+2.3 

< 19.8 

<40.7 

<36.1 

<34.7 

3C 238 

65+10 

77+12 

84+12 

<92 

<283 

266+40 

<4.0 

<8.2 

< 13.4 

< 14.2 

< 17.3 

3C 239 

96+14 

111+17 

130+20 

142+21 

651+98 

1450+218 

<7.3 

< 15.3 

< 19.1 

< 14.9 

<21.0 

3C 241 

92+14 

101+15 

116+17 

161+24 

389+58 

591+89 

7.9+1.1 

<8.2 

< 15.5 

< 12.5 

< 17.6 

3C 245 

1420+213 

1900+285 

3350+503 

5270+790 

10400+1560 

20400+3060 

47.5+2.4 

35.3+6.4 

35.2+5.9 

< 16.1 

<22.3 

3C 249 

54+8 

52+8 

42+6 

47+7 

194+29 

<516 

<3.2 

< 10.7 

< 11.5 

< 10.4 

< 14.4 

3C 250 

61+9 

59+9 

46+7 

29+4 

162+24 

< 147 

<3.7 

< 10.7 

< 11.8 

< 15.0 

< 13.9 

3C 252 

225+34 

382+57 

787+118 

1390+209 

3900+585 

7000+1050 

21.4+2.4 

<21.3 

< 17.2 

< 17.4 

<23.1 
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Table 2. continued. 


Object 

(1) 

F3.6//111 

OtJy) 

(2) 

F4.5Aim 

OJy) 

(3) 

F.5.8/7111 

OJy) 

(4) 

Fs.0/7M1 

GuJy) 

(5) 

Fl6/jm 

(m Jy) 

(6) 

F24/JII1 

C/tJy) 

(7) 

F?0/jm 

(mJy) 

(8) 

F1 60/vm 

Only) 

(9) 

F250/im 

(mJy) 

(10) 

F350//H1 

(mJy) 

(11) 

FsOfyim 

(mJy) 

(12) 

3C 255 

85+13 

86+13 

57+9 

22+3 

< 116 

<241 

<3.2 

<8.4 

< 14.5 

< 17.0 

< 19.1 

3C 256 

34+5 

37+6 

43+7 

75+11 

743+111 

1900+285 

17.8+2.3 

31.9+5.3 

28.2+6.9 

< 19.1 

<21.6 

3C 257 

85+13 

111+17 

194+29 

322+48 


1360+204 

8.1 + 1.0 

15.6+2.5 

33.1+4.7 

31.8+5.5 

32.3+8.6 

3C 266 

68+10 

73+11 

45+7 

102+15 

370+56 

980+147 

7.6+2.4 

29.4+4.1 

19.5+5.6 

< 15.9 

< 18.9 

3C 267 

153+23 

218+33 

414+62 

739+111 

2370+356 

3730+560 

15.3+2.3 

< 12.7 

< 16.8 

< 14.1 

< 19.1 

3C 268.4 

1060+159 

1560+234 

2220+333 

3330+500 

7580+1137 

11600+1740 

30.3+2.1 

< 17.0 

< 18.7 

< 18.2 

<24.3 

3C 270.1 

606+91 

944+142 

1430+214 

2260+339 

3910+587 

5470+821 

30.0+2.2 

38.9+4.7 

27.9+5.5 

< 17.5 

< 19.6 

3C 280.1 

378+57 

512+77 

777+116 

1170+176 

1680+252 

2160+324 

<5.0 

<9.4 

< 17.0 

< 13.4 

< 17.2 

3C 294 

<93 

< 103 

68+10 

67+10 


348+52 

<6.6 

<22.1 

< 16.6 

< 14.5 

<21.9 

3C 297 

119+18 

126+19 

122+18 

121+18 

<288 

432+65 

12.6+1.2 

15.4+2.4 

24.5+4.3 

< 13.8 

< 17.2 

3C 298 

1600+240 

2390+359 

3710+556 

5510+827 

9160+1374 

12600+1890 

78.7+2.4 

96.8+4.3 

96.0+6.9 

51.5+5.9 

24.1+6.7 

3C 305.1 

181+27 

282+42 

495+74 

972+146 

2410+362 

2490+374 

24.0+2.3 

40.4+4.3 

34.9+6.0 

< 18.1 

< 18.9 

3C318 

343+51 

427+64 

571+86 

806+121 

1960+294 

3400+510 

18.3+2.3 

43.8+5.9 

42.9+6.1 

b 

b 

3C 322 

128+19 

135+20 

94+14 

120+18 

411+62 

804+121 

<5.0 

< 11.0 

< 19.0 

< 17.5 

< 21.9 

3C 324 

165+25 

160+24 

178+27 

450+68 

2580+387 

2820+423 

23.5+2.3 

31.7+5.6 

21.0+6.2 

< 18.2 

< 21.4 

3C 325 

472+71 

565+85 

708+106 

1200+180 

1990+299 

3030+455 

<7.6 

< 18.4 

< 15.4 

< 16.1 

< 18.3 

3C 326.1 

29+4 

34+5 

26+4 

72+11 

829+124 

1430+215 

<6.7 

< 14.5 

29.8+5.9 

21.0+6.9 

<20.8 

3C 356 

108+16 

110+16 

122+18 

434+65 

2270+341 

4060+609 

11.6+2.5 

19.7+4.9 

< 18.4 

< 15.3 

<20.4 

3C 368 

126+19 

112+17 

112+17 

210+32 

1370+206 

3250+488 

29.9+2.0 

61.5+4.8 

44.4+7.4 

23.8+6.2 

<21.3 

3C418 

1130+170 

1630+245 

2470+371 

3900+585 

6680+1002 

13600+2040 

95.1+2.6 

200.1 + 16.1 

173.8+40.0 

259.2+43.7 

387.4+28.6 

3C 432 

420+63 

526+79 

857+129 

1490+224 

2710+407 

3940+591 

<7.3 

< 13.8 

34.7+5.1 

25.2+5.1 

30.7+5.8 

3C 437 

82+12 

85+13 

97+15 

80+12 

384+58 

941+141 

<6.4 

< 18.4 

< 17.0 

< 14.3 

< 18.9 

3C 454.0 

339+51 

481+72 

811+122 

1220+183 

2490+374 

4150+623 

15.7+2.2 

39.5+5.2 

31.0+5.6 

35.3+5.2 

28.0+7.2 

3C 454.1 

77+12 

76+11 

112+17 

135+20 

612+92 

1500+225 

13.7+2.5 

37.0+4.7 

50.2+8.7 

26.7+10.4° 

<50.0 

3C 469.1 

160+24 

244+37 

509+76 

1090+164 

3270+491 

1970+296 

10.9+2.3 

24.3+4.6 

<20.9 

< 19.1 

<21.8 

3C 470 

50+7 

75+11 

72+11 

266+40 

1510+227 

2650+398 

16.0+2.7 

29.3+5.1 

48.0+6.5 

36.3+5.2 

<21.5 

4C 13.66 

24+4 

24+4 

21+3 

18+3 

<260 

276+41 

<5.5 

< 13.4 

< 15.6 

< 13.9 

< 17.9 

4C 16.49 

329+49 

420+63 

573+86 

743+111 

1070+161 

1830+275 

<5.0 

< 17.3 

< 18.5 

< 18.7 

<23.2 


Notes. 

(a) Less than 3<r detection entering our SED fitting routine. 

<w Photometric measurement hindered by the presence of a nearby source. 
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Table 4. Physical parameters estimated from the SED fitting. 


Object 

Lagn 

Lsf 

SFR 

M d 

Td 


(10 12 l g ) 

(10 i 2 L o ) 

(10 2 M q yr -1 ) 

(10 8 Me) 

(K) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Objects detected in at least three Herschel bands 

3C 002 

2-lti 

o.618:t 

i-i28:l 

i.i28:? 

31.62i? 

3C014 


1 .318:2 

2-328:1 

0.92oi 

37.52?;? 

3C 068.2 

6-628:1 

2.128:! 

3-728:1 

7-021:8 

28.52?;? 

3C 119° 

4-4^8:I 

0.928:? 

1-528:1 

0-728:1 

37-62?;? 

3C 124 

3.7-8:! 

1-528:1 

2-728:1 

2-628:8 

31.92U 

3C 190 

9.o-8i 

2-728:1 

4-728:? 

4-928:1 

31.7-8J 

3C 205 

oooo 
do 
+ 1 

O ) 

(N 

4-128:1 

7-028:1 

1-828:1 

40.721:8 

3C 222 

1-518:1 

2-528:1 

4-328:1 

5-428:1 

30.628:8 

3C 245 

11.218:1 

0-828:1 

1-428:1 

2-428:8 

90 Q+3.4 
420.5_2o 

3C 256 

6-918:1 

2-628:1 

4-528:1 

0-928:1 

43.0211 

3C 257° d 

8-018:6 

5-428:1 

9-228:1 

3-428:1 

38.228:8 

3C 266 


1-528:1 

2-728:1 

0-428:1 

44.22i9 

3C 270.1 

i2.7i8i 

2.328:1 

3-928:5 

0-828:1 

43.52!i 

3C 297 

9 9+0.1 

z.z_ 02 

0.928:1 

1-628:1 

6.52?:? 

24.921:8 

3C 298* 

29.218J 

5-428:1 

9.328:1 

9 c+0.3 
J.o_ 04 

37.5281 

3C 305.1“ 

3-428:1 

1-328:? 

9 9+0.4 
z>z -0.2 

1-428:1 

34.7218 

3C318 

7-628:? 

3.428:1 

5-828:1 

1-728;! 

39-62?;? 

3C 324 

3.618:1 

1-028:1 

1-828:1 

0-528:1 

39.6-55 

3C 368* 

3.318:1 

2-028:1 

3-528:1 

1-428:1 

37.32?j 

3C 432* 

9.318:2 

2-428:1 

4-228:1 

6.42!i 

29.821:8 

3C 454.0 

10.5185 

3-628:1 

6-228:1 

2-328:8 

38.32?;8 

3C 454.1 

5.218:? 

4-428:1 

7-528:1 

2-428:1 

39.32?;8 

3C 470* 

5-828:1 

3-028:1 

5-128:1 

4.528;? 

32.821:8 

Objects detected in fewer than three Herschel bands 

3C 009 

17.1 

< 1.8 

< 3.1 



3C013 

4.0 

< 1.2 

<2.0 



3C 036“* 

< 1.2 

<0.5 

<0.9 



3C043 

2.8 

< 1.2 

<2.0 



3C065“* 

< 1.4 

<0.6 

< 1.0 



3C 068.1“ 

8.6 

<0.7 

< 1.1 



3C 173* 

0.4 

<0.3 

<0.5 



3C 181“ 

5.2 

<0.8 

< 1.5 



3C 186“ 

4.7 

<0.5 

<0.8 



3C 191 

17.0 

< 1.7 

<3.0 



3C 194* 

<0.7 

<0.6 

<0.9 



3C 204“ 

<4.8 

<0.5 

<0.9 



3C 208.0 

< 3.6 

<0.6 

< 1.0 



3C 208.1“ 

< 1.6 

<0.4 

<0.7 



30 210“ 

5.2 

< 1.9 

<3.3 



3C212 

6.2 

<0.8 

< 1.3 



3C 220.2 

6.3 

<0.7 

< 1.2 



3C 225A* 

< 1.5 

< 1.8 

<3.1 



3C 230 

3.1 

< 1.9 

< 3.3 



3C 238* 

<0.7 

<0.5 

<0.9 



3C 239* 

<3.0 

< 1.7 

<2.9 



3C 241* 

2.1 

<0.8 

< 1.4 



3C 249* 

<0.9 

<0.8 

< 1.5 
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Table 4. continued. 


Object 

Lagn 

Lsf 

SFR 

M d 

T d 


(10 12 l g ) 

(1 O i2 L 0 ) 

(10 2 M Q yr 1 ) 

(10 8 Mo) 

(K) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

3C 250 

<0.6 

<0.5 

<0.9 



3C 252“ 

4.5 

<0.6 

< 1.0 



3C 255* 

<0.5 

<0.5 

<0.9 



3C 267“ 

3.1 

<0.4 

<0.7 



3C 268.4 

16.2 

<0.8 

< 1.4 



3C 280.1“ 

<5.4 

< 1.1 

< 1.9 



3C 294* 

<2.3 

< 1.7 

<2.9 



3C 322 

< 1.5 

< 1.3 

<2.2 



3C 325 

<2.5 

<0.7 

< 1.2 



3C 326.1 

<3.7 

<3.0 

<5.2 



3C 356 

1.9 

<0.8 

< 1.3 



3C 437* 

< 1.6 

< 1.2 

<2.1 



3C 469.1“ 

3.3 

< 1.3 

<2.2 



4C 13.66 

< 1.0 

<0.8 

< 1.4 



4C 16.49“ 

<2.4 

<0.9 

< 1.5 




Notes. (1) Name of object; (2) IR luminosity (integrated between 1 /rm and 1000 pm) of the AGN powered dust emission, i.e. torus component and 
sum of torus and hot dust components for RGs and QSRs, respectively; (3) IR luminosity (integrated between 8 pm and 1000 pm) of the modified 
blackbody component (fi = 1.6) presumably powered by star formation activity in the AGN host galaxy; (4) star formation rate determined from 
the IR luminosity in (3), using the calibration derived in lKennicuttl d 19981) ; (5) mass of the modified blackbody component; (6) temperature of the 
modified blackbody component. 

(a) p or this rac po galaxy, a blackbody component (1300 K) was included in the SED fitting. 

(W This object was included in the stacking of non-detected radio galaxies. 

(c ' This object was included in the stacking of non-detected quasars. 

The emissivity index of the cold dust component, fl, was estimated for this object. 
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Appendix A: UV/visible SEDs of quasars 

The emission from powerful quasars at UV/visible wavelengths 
comes predominantly from the hot accretion disk. In quasar SED 
fitting, this component is typically represented by a power-law 
(PL), extended to the NIR by introducing a PL index of -2 
(F v oc v 2 ) (e.g. lLeipski et al.ll2013l) In the current work, we do 
not fit the rest-frame UV/visible SEDs of quasars because it has a 
negligible influence on the LIR part of the SED. We demonstrate 
this using two objects that have good coverage in the UV/visible: 
3C 454.0 and 3C 002. 

Firstly, we consider 3C 454.0 (upper panels of Fig. ED , 
for which we have SDSS photometry ( ugriz ) in addition to the 
Spitzer and Herschel photometry, and compare the results when 
we include/exclude the PL component. In a first attempt, we fit 
the data using a PL component in addition to the 1300 K black- 
body, torus, and star-formation-heated cold dust components. 
Given the redshift of 3C 454.0, strong Mgll emission at rest- 
frame ~ 0.3 pm, typical for optical spectra of quasars, signifi¬ 
cantly contributes to its broad-band i and z photometry. The SED 
fitting including the PL component was therefore performed us¬ 
ing only the u, g, and r SDSS photometry. In a second attempt, 
we fit only data at rest-frame wavelengths longer than 2 pm, ex¬ 
cluding the PL component from the fit. As seen in Fig. IA.ll the 
best-fit SEDs longward of 2 pm rest-frame wavelength are very 
similar for the two fitting approaches. The only difference is the 
flux normalization of the hot dust component, which turns out to 
be slightly lower when including the UV/visible part of the SED. 
Understandably, this is a result of the PL component contribut¬ 
ing to the emission in the shorter Spitzer wavelengths. While the 
fitting procedures prefer different torus models, the luminosity of 
this component changes by less than 10%. Most importantly, the 
most relevant physical parameters for the current work: the star 
formation rate, and the temperature and mass of the cold dust 
component remain within 10%. 

Secondly, we consider 3C 002 (lower panels of Fig. IA. 1 1 . 
for which in addition to SDSS, Spitzer and Herschel photome¬ 
try, we have photometry from 2MASS. We fit the full SED of 
this object, after discarding data points potentially contaminated 
by strong emission lines. Once again we reach the same results 
as in the case for 3C 454.0, thus we conclude that a detailed 
treatment of the emission in the UV/visible part of the quasars’ 
SEDs has a negligible effect on the results inferred from the 
IR part of the SED. Furthermore, many 3CR quasars lack good 
UV/visible/NIR photometric data, making the inclusion of the 
PL component to the fitting procedure impossible. Consequently, 
the results presented in our work were obtained using the best-fit 
SEDs without fitting a PL to the UV/visible part of the SEDs of 
quasars. 


data. 

3C 208.1 - This source requires an additional hot dust compo¬ 
nent to better fit the data. 

3C 210 - Bright nearby object present in SPIRE bands. De¬ 
blending is not possible. We report no measurements in the 
SPIRE bands. Photometry in PACS 160 pm was performed with 
an aperture of 6" radius. This source requires an additional hot 
dust component to better fit the data. IRAC 8 pm and IRS 16 pm 
data are not well fitted. 

3C 222 - IRS 16 pm data point is removed from the fitting 
procedure. 

3C 230 - The Spitzer photometry probably includes a star 
located very close to the radio galaxy. IRAC 3.6 pm, IRAC 
4.5 pm, and IRAC 5.8 pm data points are removed from the 
fitting procedure. 

3C 252 - This source requires an additional hot dust component 
to better fit the data. 

3C 255 - The IRAC points are not well-fitted with the blackbody 
component representing emission from old stars in the host 
galaxy. 

3C 257 - This source requires an additional hot dust component 
to better fit the data. 

3C 267 - This source requires an additional hot dust component 
to better fit the data. 

3C 305.1 - This source requires an additional hot dust compo¬ 
nent to better fit the data. IRAC 8 pm and IRS 16 pm data are 
not well fitted. 

3C 318 - Bright nearby object present in PACS and SPIRE 
bands. Deblending is not possible for SPIRE 350 pm and SPIRE 
500 pm. We report no measurement in these two bands. 

3C 324 - IRAC 8 pm and IRS 16 pm data are not well fitted. 

3C 418 - The only flat-spectrum-core-dominated object within 
the high-z 3CR sample; completely synchrotron dominated. We 
do not include this source in the analysis. 

3C 454.1 - Diffuse emission present in SPIRE maps. 

3C 469.1 - The PACS 160 pm flux density might be slightly 
contaminated by a nearby object. This source requires an 
additional hot dust component to better fit the data. IRAC 8 pm 
and IRS 16 pm data are not well fitted. 


Appendix B: Comments on individual objects 

3C 036 - This source requires an additional hot dust component 
to better fit the data. 

3C 043 - There is a bright nearby object dominating the emission 
in the PACS 160 pm band. We report no measurement in this 
band. 

3C 065 - This source requires an additional hot dust component 
to better fit the data. 

3C 068.2 - IRS 16 pm data point is removed from the fitting 
procedure. 

3C 119 - Diffuse emission present in the PACS 160 pm and 
SPIRE maps owing to the object’s low galactic latitude. This 
source requires an additional hot dust component to better fit the 
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Fig. A.l. Spectral energy distributions of two quasars with good visible/UV detections. Left panels include a power-law component (dash-dotted 
pink) representing the emission from the hot accretion disk. Right panels exclude that component, and show the components considered when 
fitting the quasars in this work. Other individual components as described in Fig. [2] 
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Appendix C: Best-fit SEDs of objects detected in at 
least three Herschel bands 

Appendix D: Best-fit SEDs of objects detected in 
fewer than three Herschel bands 

Appendix E: Postage stamps 
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Fig. C.l. Spectral energy distributions of the 24 objects detected in at least three Herschel bands. Individual components as described in Fig. [2] 
3C 418 was not included in the analysis, as discussed in Sect. 13.31 
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Fig. D.l. Spectral energy distributions of the objects detected in fewer than three Herschel bands. Individual components as described in Fig. [2] 
Red circles denote 3cr upper limits taken to be tentative detections when calculating upper limits of physical parameters, as explained in Sect. [3] 
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Fig. D.l. Continued. 




























































































































Fig. E.l. Postage stamps of the high-z 3CR objects studied in this work. From left to right: IRAC 3.6 /im, IRAC 4.5 /im, IRAC 5.8 pm, IRAC 
8 /mi. IRS 16 /mi, MIPS 24 /mi. PACS 70 /mi, PACS 160 /an, SPIRE 250 /mi, SPIRE 350 /mi, and SPIRE 500 //in bands, respectively. Each image 
shown here has dimensions of 2x2 arcmin. The circle (10" in radius) indicates the central position of the source. 
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Fig. E.l. Continued. 
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